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We present new data on Cherenkov light observations obtained during the period 1994-2009, after a modern- 
ization of the Yakutsk EAS array. A complex analysis of Xmax and its fluctuations <r(xmax) was performed 
over a wide energy range. With the new data, according to QGSJet II model, an estimation was made of the 
cosmic ray mass composition for Eo ~ 10 17 — 3 X 10 19 eV. The result points towards a mixed composition 
with a large portion of heavy nuclei at Eo ~ 10 17 eV and the dominance of light nuclei at Eo ~ 10 19 eV. The 
analysis of a(x max ) energy dependence for the same energies qualitatively confirms this result. The shape of 
the x max distribution at fixed energy 10 18 eV is analysed to make more precise conclusions on cosmic ray mass 
composition. 



1. Introduction 

The Yakutsk EAS array effectively covers the en- 
ergy domain from 10 15 eV to 5 x 10 20 eV by measuring 
charged particles, muons with e t h r > 1 GeV- sec 9 and 
Cherenkov light emission. In the same energy region 
two irregularities of cosmic ray (CR) energy spectrum 
are observed — the knee (3 x 10 15 eV) and the an- 
kle (8 x 10 18 eV). The nature of these irregularities is 
yet to be explained. From the work by Berezinsky et 
al [l[ it follows that such a spectrum shape could be 
a consequence of a transition from galactic to extra- 
galactic components in the total CR flux, i.e. there 
must be some region in the CR spectrum where the 
intensities of these fluxes become equal to each other 
and then decrease. The boundaries of such a transi- 
tion region is yet unknown and represents a subject 
of research at many large EAS arrays. The physics 
of this phenomenon is tightly connected with particle 
drift in the magnetic fields of our Galaxy and outer 
space. Since nuclei of different masses behave differ- 
ently in a magnetic field, then the CR composition in 
different points of space will differ. Therefore, it is 
possible to solve the problem of the transitional re- 
gion by measuring the CR mass composition in EAS 
related experiments. 

It is a known fact that the depth of shower maxi- 
mum (x max ) and fluctuations in EAS development are 
sensitive to the atomic number of the primary par- 
ticle and for this reason they are used to estimate 
the CR mass composition 0t3|- It is especially im- 
portant for the ultra-high energy region where direct 
measurements of mass composition are impractica- 
ble. For example, works by Dyakonov et al [!, H, 0] 
utilized single characteristics and their combinations: 
(xmax), o-(x max ) and dx rn ax/dlgE . These works pro- 
vided initial estimations of inelastic interaction cross- 
sections at ultra-high energies. 

In this paper we present the data on longitudi- 
nal EAS development reconstructed from Cherenkov 
emission data. These data were obtained after mod- 
ernization of the Yakutsk array when the accuracy of 



the main EAS characteristics increased compared to 
the previous series of observations. A precise knowl- 
edge of the mass composition together with the energy 
spectrum plays a major role in understanding CR as- 
trophysics 8]. In this sense, engaging the maximal 
possible number of composition-sensitive EAS charac- 
teristics increases the reliability of CR chemical com- 
position estimation. It is important to consider not 
only the mean EAS parameters, e.g x max , muon con- 
tent p c h but also their fluctuations in given energy 
intervals 0, In order to minimize the latter, it is 
also a good idea to analyze them at fixed energies. 



2. Technical aspects of longitudinal EAS 
development characteristics estimation 

The determination of x max in individual showers 
is based on methods developed at the Yakutsk ar- 
ray and utilize the measurements of EAS Cherenkov 
light emission at different core distances. In the 
first method, x max is determined by the parameter 
p = lg Q200 / Q550 (a relation of Cherenkov light fluxes 
at 200 and 550 m from the core); the second — in- 
volves the reconstruction of the EAS development cas- 
cade curve, using the Cherenkov light lateral distribu- 
tion function and a reverse solving the third is 
based on half- width and half-height of Cherenkov light 
pulses recorded at 200 m from the core; the fourth 
method includes recording the Cherenkov track with 
several detectors based on camera-obscura located at 
300 — 500 m from the array center [l3| . 

Various factors affect the methods mentioned 
above: the way the showers are selected, precision of 
core location, atmosphere transparency, mathemati- 
cal methods used to calculate parameters of approx- 
imated functions, hardware-related fluctuations and 
so on. The influence from single and composite fac- 
tors on physical results of the Yakutsk array opera- 
tion was calculated either with full simulation of the 
measurement procedure or empirically estimated dur- 
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Figure 1: Individual EAS events detected in the Yakutsk 
experiment during 1994 - 2009. 



ing special methodical experiments. For instance, an 
estimation of hardware-related errors was performed 
by the analysis of two nearby detectors that measure 
charged particles, muons and Cherenkov light emis- 
sion [(| [lj| . The accuracy of x max determination in 
individual showers was estimated in simulating EAS 
characteristics measurements at the array involving 
Monte-Carlo methods and amounted to 30 — 45 g/cm 2 , 
35 — 55 g/cm 2 , 15 — 25 g/cm 2 , 35 — 55 g/cm 2 respec- 
tively for the first, second, third and fourth meth- 
ods. The total error of x max estimation included er- 
rors associated with core location, atmospheric trans- 
parency during the observational period, hardware 
fluctuations and mathematical methods used to cal- 
culate main parameters. 



3. Mean depth of maximum shower 
development 

Figured] demonstrates a cloud of points in the x max 
distribution for showers with energy above 10 17 eV. 
These data were obtained using all four methods and 
reflect an alteration of x max towards lower atmosphere 
depths with increasing energy. Figure [5] shows x max 
values averaged over energy intervals together with 
the data from other experiments. On the same pic- 
ture results of different hadron models calculations are 
shown. All experimental data coincide within experi- 
mental errors and demonstrate an irregular shift with 
energy. Up to 3 x 10 18 eV E.R. has value 60-80 g/cm 2 
and within the interval of 3 x 10 18 — 5 x 10 19 eV it 
equals 40 — 60 g/cm 2 . This might be interpreted as a 
possible alteration in mass composition at very high 
energies. A comparison with calculations reveals a 
tendency for a light nuclei abundance starting from 
5 x 10 17 eV to 2 x 10 19 eV and some abundance above 
2 x 10 19 eV. 



Figure 2: Energy dependence of i mal . Filled circles 
represent the Yakutsk data, open circles — CASA-MIA, 
squares — AUGER data, blue triangles — preliminary 
results of the Telescope Array experiment [lg| . Solid 
lines — results obtained with QGSJet II, dashed — 
EPOS 1.6, point line — SIBYLL 1.62 



4. Fluctuations of x, 



Fluctuations of x max play a huge role in EAS lon- 
gitudinal development as they are associated with the 
point of first interaction (and, hence, with the inelas- 
tic interaction cross-section, <7A-air), energy transfer to 
secondary hadron particles (inelastic coefficient -Ki ne i) 
and, to a great extent, depend on the kind of primary 
particle initiating a shower. So, the amount of fluc- 
tuations measured in different energy intervals could 
characterize the CR mass composition at a given en- 
ergy and on the whole determine the dynamics of its 
change with the energy of the primary particle. Fig- 
ure [3] demonstrates the energy dependence of o~(x max ) 
obtained at the Yakutsk array and for comparison the 
same figure shows the HiRes data [ItJ • The data from 
the HiRes experiment virtually reproduces the data 
from Yakutsk but have a slight tendency of o~(x max ) 
change: a small increase in the region of 10 17 — 10 18 eV 
and a decrease at 2 x 10 18 — 5 x 10 19 eV. The curves rep- 
resenting simulation results, obtained with QGSJetOl, 
QGSJet 11 and SIBYLL models, are also shown on this 
figure. Calculations were performed for proton, he- 
lium nuclei, CNO group and iron nuclei. Comparison 
with experimental data has shown that the CR com- 
position in this energy region is mixed with a dom- 
inance of protons and helium nuclei. It should be 
pointed out that according to Figure [3l the portion of 
heavy nuclei in the CR flux energy above 2 x 10 18 eV 
is small and helium and CNO-group nuclei might play 
a significant role. We came to the same conclusion [J] 
where the shape of the x max distribution was analyzed 
within the framework of the QGSJetOl model at fixed 
energies 10 18 eV and ~ 10 19 eV (see Figure S]). 
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Figure 3: Fluctuations of the depth of maximum EAS 
development: filled squares — Yakutsk data, open 
squares — HiRes data, open triangles — data from 
Pierre Auger Observatory. Straight line — results 
obtained with QGSJetOl, dashed line — QGSJet II, 
dotted line — SIBYLL 1.62 for various primary nuclei 
(see Abbasi et al [r^l 



Figure 5: Mean mass number of primary particle as a 
function of energy. Circles represent the Yakutsk data, 
triangles — HiRes data, squares — results obtained at 
Auger observatory, blue empty triangles — preliminary 
data from the Telescope Array experiment 11611 . dotted 
line — computational results by Berezhko 
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Figure 4: x ma x distribution at fixed energy 10 18 eV. 
Solid line represents the Yakutsk data 
(8 x 10 17 < E < 2 X 10 18 eV, (E ) = 1.0 x 10 18 eV, 
857 events); dotted line — QGSJetOl for mixed 
composition (70% p, 30% Fe); dashed line — QGSJetOl 
for primary protons, solid grey line — QGSJetOl for 
CNO group nuclei, dash-dotted line — QGSJetOl for 
iron nuclei (see Knurenko et al Q) 

5. Cosmic ray mass composition 

Figure [5] displays the mean natural logarithm of the 
CR atomic number (In A) concluded from the x max 
data from four experiments — Yakutsk, HiRes, Auger 
and Telescope Array [l6| . For (In A) derivation x max 
values were utilized, obtained in simulations within 
the framework of the QGSJet II models for proton and 
iron nuclei. The (In A) value was calculated according 
to the relation proposed by Horandel fl8j : 

{lnA) = lZ aX ~ jr -1*56 (1) 



At first glance, all data reveal a tendency to change 
(In A) with energy. For instance, in the energy interval 

2 x 10 17 - 3 x 10 18 eV, the value of (In ,4} drops from 

3 to 1.3 and above 10 18 eV a slight increase is noted. 
Such a behaviour is close to the "dip" -scenario from 
the work by Berezhko fl9| , where two peaks are seen 
in the energy dependence of (In A). The first one, at 
~ 10 17 eV, corresponds to the ending of the galactic 
component, the second — at 10 19 — to the start of 
CR intensity change due to GZK-cutoff. 

However, there is still a significant data dispersion 
in this energy region due to poor event statistics. 
Thus, the reliability of our statement is quite limited. 
For a more precise conclusion on ultra-high energy 
cosmic rays origin, a few conditions must be fulfilled: 
improved statistics, improvement of x max estimation 
precision, adaptation of a single hadron interaction 
model that well describes experimental data and in- 
volving several alternative methods for x max evalua- 
tion. 



6. Conclusions 

Thus, according to all the data reviewed above, 
within the framework of the QGSJet hadron inter- 
action model it is reasonable to speculate that the 
primary cosmic ray mass composition alters during 
the energy transition from 10 17 eV to 5 x 10 18 eV. 
At E > 5 x 10 18 eV ~ 70% of cosmic rays consist 
of protons and helium nuclei. The content of other 
nuclei in the region of the ankle of the spectrum does 
not exceed ~ 30 %. A large portion of protons and he- 
lium nuclei in the primary CR near the ankle is most 
likely associated with a significant contribution from 
particles arriving from outside our Galaxy. In such 
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Table I New data from the Yakutsk array 



E , eV 


1.5 x 10 17 


2.0 x 10 17 


3.0 x 10 17 


4.5 x 10 17 


5.8 x 10 17 


8.6 x 10 17 


1.2 x 10 18 


1.7 x 10 18 


2.4 x 10 18 


%max 


632 


648 


655 


670 


687 


685 


700 


719 


723 




5 


5 


6 


6 


7 


7 


7 


7 


8 


E , eV 


3.5 x 10 18 


5.0 x 10 18 


8.0 x 10 18 


1.1 x 10 19 


1.6 x 10 19 


2.4 x 10 19 


3.5 x 10 19 


5.1 x 10 19 




%max 


738 


765 


761 


786 


780 


778 


785 


808 




CT^Xmax) 


8 


9 


11 


13 


18 


21 


26 


25 





a case the region of transition from galactic to extra- 
galactic component might be in the energy interval 
10 17 — 10 19 eV. The problem of mass composition al- 
tering above 10 19 eV remains unresolved due to poor 
event statistics. 
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